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Abstract Mayall II = Gl is one of the most luminous globular clusters (GCs) in M31. 
Here, we determine its age and mass by comparing multicolor photometry with theoreti- 
cal stellar population synthesis models. Based on far- and near-ultraviolet GALEX pho- 
tometry, broad-band U BVRI, and infrared JHKs 2MASS data, we construct the most 
extensive spectral energy distribution of Gl to date, spanning the wavelength range from 
1538 to 20,000A. A quantitative comparison with a variety of simple stellar population 
(SSP) models yields a mean age that is consistent with Gl being among the oldest build- 
ing blocks of M31 and having formed within ^1.7 Gyr after the Big Bang. Irrespective 
of the SSP model or stellar initial mass function adopted, the resulting mass estimates 
(of order lO^i\/0) indicate that Gl is one of the most massive GCs in the Local Group. 
However, we speculate that the cluster's exceptionally high mass suggests that it may not 
be a genuine GC. We also derive that Gl may contain, on average, (1.65 ±0.63) x IO'^Lq 
far-ultraviolet-bright, hot, extreme horizontal-branch stars, depending on the SSP model 
adopted. On a generic level, we demonstrate that extensive multi-passband photometry 
coupled with SSP analysis enables one to obtain age estimates for old SSPs to a similar 
accuracy as from integrated spectroscopy or resolved stellar photometry, provided that 
some of the free parameters can be constrained independently. 
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1 INTRODUCTION 

Globular clusters (GCs) are among the oldest bound stellar systems in the Universe, and they thus 
provide a fossil record of the earliest stages of galaxy formation and evolution. GCs are internally 
homogeneous in age and metallicity, i.e. they are simple stellar systems composed of coeval stellar 
populations. In addition, GCs are the oldest systems in our own and other galaxies for which we can es- 
timate reasonably reliable ages (and realistic uncertainties); they can thus independently provide vitally 
important information regarding the minimum age of the Universe. In a detailed study of 17 Galactic 



* E-mail: ma jun@vega . bac . pku . edu . cn 



2 



Ma et al. 



GCs, IChaboyer et al. I d 19981) used the improved Hipparcos parallaxes having just become available at 
that time to determine updated distances, and hence improved ages, of their GC sample. They con- 
cluded that the mean age of the oldest GCs i s 11.5 zt 1.3 Gyr, alth ough their age histogram (their fig. 2) 
shows a tail toward ages as old as ^ 16 Gvr. iGratton et al. 1 (l2003b obtain improved ages (and distances) 
for three Galactic GCs, NGC 6397, NGC 6752, and 47 Tuc, and conclude that the age of the oldest 
GCs is 13.4 ± 0.8 (random) ±0.6 (systematic) Gyr, in good agreement with the 3-year results from the 
Wilkinson Microwave Anisotropy Probe (WMAP). This led them to suggest that the oldest Galactic GCs 
formed within the first 1.7 Gyr after the Big Bang, at the la confidence level. We note that this is still 
fully compatible with the 5-year WMAP results constraining the age of the Universe to 13.73 ± 0.12 
Gyr (Hinshaw et al. 2008). While the ages of the oldest GCs in the Galaxy are now reasonably well 
determined, this is certainly not the case even for our nearest large neighbour, the Andromeda galaxy 
(M31). 

The most direct method for determining the age of a star cluster is by means of individual stellar 
photometry, since the main-sequence turn-off location is mostly affected by age (see, e.g., Puzia et al. 
2002b, and references therein). However, this method has only been applied to Galactic GCs and to 
GCs associated with the Milky Way's satellites (e.g., Riich 2001), in which individual stars can be re- 
solved and their photometry determined to satisfactory accuracy, to a few magnitudes fainter than the 
main-sequence turn-offs. This is difficult, if not impossible, to achieve even for GCs as close as those 
associated with M31 (but see Brown et al. 2004), at a distance of D = 772 ± 44 kpc (e.g., Ribas 
et al. 2 005). Fortunately, starting from the pioneering work of Tinsley (1968, 1972) and ISearle et al. I 
dl973h . evolutionary population synthesis modeling has become a powerful tool to address many out- 
standing problems in astrophysics, from determining the ages of star clusters to interpreting integrated 
spectrophotometric observations of galaxies. Therefore, extragalactic GC ages can, in general, also be 
inferred from composite colors and/or integrated spectroscopy. 

The evolution of GCs is usually modeled by means of the simple stellar population (SSP) approx- 
imation. An SSP is a single generation of coeval stars formed from the same progenitor molecular 
cloud (thus implying a single metallicity), and governed by a given stellar initial mass function (IMF). 
GCs are ideal templates to test the compatibility between the population synthesis models and reality. 
iBarmbv & Huchral ( l2000l) compared the predicted SSP colors of three stellar population synthesis mod- 
els to the intrinsic broad-band U BVRIJHK colors of Galactic and M31 GCs, and concluded that the 
best-fitting models match the cluster colors very well. Subsequently, many authors have use d SSP mod- 
eling to determine the parameters of cluster populations. For instance. lde Griis et al. I (l2003al) determined 
the ages and masses of star clusters in the fossil starburst region B of M82 by comparing the observed 
cluster spectral energy distribut ions (SEDs) with both the STARBURST99 SSP models (Leitherer et al. 
1999) and those developed bv iBruzual & Chariot I ( |2000|) , based on Hubble Space Telescope (HST) 
observations fro m th e blue optical t o the near-infrared (NIR) (see also de Grijs et al. 2003b, 2003c); 
iBik et al. I ( l2003b and lBastian et~an ( l2005b estimated ages, initial masses and extinction values for M5 1 
star cluster candidates by comparing the STARBURST99 and the Frascati models (Romaniello 1998) for 
instantaneous star formation with t he observed SED s bas ed on HST/WFPC 2 observations in six broad- 
band and two narrow-band filters. iMa et al.1 (|2006|) and iFan et al. I (l2006l) obtained age estimates for 
M3 1 GCs by fitting theoretical stellar population synthesis models (Bruzual & Chariot 2003, henceforth 
BC03) to their photometric measurements in a large number of inter mediate- and broad -band passbands 
from the optical to the NIR. Based on the same method and models. iMa et al. I (12007 ah constrained the 
age of the M31 GC S3 12, using multicolor photometry from the near-ultraviolet (NUV) to the NIR, 
to 9.5lQ gg Gyr. S312 is among the first extragalactic GCs for which the age was estimated accurately 
using main-sequence photometry, i.e.. lBrown et al. I (12004') estimated the age of S312 at lOl^i'^ 
means of a quantitative comparison with the isochrones of VandenBerg et al. (2006J. This was based on 
their analysis of the cluster's color-magnitude diagram (CMD) below the main-sequence turn-off using 
extremely deep images obtained with the HST/Advanced Camera for Surveys (ACS). 

It is a common misconceptio n that spectroscop i c age estimates are always much better than those 
based on broad-band photometry. ISchweizer et al. I (|2004|) recently showed convincingly that spectro- 
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scopic age determinations are not necessarily better o r more accurate than photometrically obtained 
ages, at least in the age range from - 100 - 500 Myr lAnders et al.1 ( |2004 published a detailed the- 
oretical investigation of the accuracy of retrieved star cluster properties, including their ages, based on 
sophisticated fits of SSP models to observed broad-band SEDs spanning varying wavelength ranges. 
They concluded that if one has access to as large a wavelength range as possible, ideally including both 
ultraviolet and NIR data points, the resulting age estimates are reasonably accurate, even for ages as old 
as ^ 10 Gyr - particularly if one or more of the other free parameters (e.g., metallicity or extinction) 
can be constrained reliably and independently (see Anders et al. 2004, e.g., their fig. 14). We will use 
this promising approach as our basic premise in this paper. 

Of the Local Group members, M31 is particularly important as it provides a direct comparison 
with our own Galaxy. In addition, it contains a large number of GCs and GC candidates, including 
496 genuine, 367 probable, and 301 possible GCs (Kim et al. 2007). The M31 GC system is among 
the extragalactic GC systems studied most often (Harris 1991; Brodie & Strader 2006). As one of the 
brightest M31 GCs, Mayall II = Gl has attracted much scientific interest (see, e.g., Barmby et al. 2007; 
Ma et al. 2007b, and references therein). 

In this paper, we first determine the age and mass of Gl by comparing observational SEDs (§2) 
with population synthesis models (§3). We will use the lessons learned from studies of broad-band 
photometric SED fits to minimize the associated uncertainties. We discuss our results along the way, as 
appropriate, and provide a summary in §4. 

2 ULTRAVIOLET, OPTICAL, AND INFRARED OBSERVATIONS OF Gl 

2.1 Historical observations 

Gl was first detected bv lMavall & Eggenl (Il953b (their No. 2, and hence referred to as Mayall II), who 
searched for nebulous objects associated wit h M3 1 using a 6° x 6° Palomar 48-inch Schmidt plate taken 
in 1 948 and centered on M3 1 . Subsequently, ISargent et al. I d 19771) rediscovered the cluster (their No. 1 , 
i.e., Gl) based on their survey of 29 plates associated with the general field of M31, which had been 
obtained at the //2.7 prime focus of the Kitt Peak National Observatory's (KPNO) 4-m telescope. The 
cluster is located in the halo of M31, at a projected distance of about 40 kpc from the galaxy's nucleus 
(see Meylan et al. 2001). 

2.2 GALEX ultraviolet, optical broad-band, and 2MASS NIR photometry 

Although the cluster is generally believed to be among the oldest GCs in M3 1 , to the best of our knowl- 
edge there is no CMD-based or spectroscopic age estimate available in the literature to date. The lack 
of a CMD-based age estimate is due to the challenges associated with probing the cluster's CMD down 
to below the main-sequence turn-off. The current deepest CMD of the cluster (Meylan et al. 2001) does 
not reach these faint levels. Although both integrated and spatially-resolved spectra of the cluster are 
available (e.g., Huchra et al. 1991; Gebhardt et al. 2005; Cohen 2006), they have thus far not been used 
to determine an age for Gl. This may be partially due to the Umited wavelength range covered by most 
of these spectra, and the difficulties one faces when trying to constrain ages in the regime beyond ~ 10 
Gyr (see below). 

To constrain the age of Gl accurately, with the smallest uncertainty allowed by the observa- 
tional data, w e use as many photometric data points covering as large a wavelength range as possible. 
iKavirai et al. I ( 12007.) showed that the combination of far (FUV) and near-ultraviolet photometry with 
optical obse rvations in the s tandard broad bands enables one to efficiently break the age-metallicity 
degeneracy; IWorthevl(ll994 showed that the age-metallicity degeneracy assoc iated with optic al broad- 
band colors is Aage/AZ ~ 3/2 (also see MacArthur et al. 2004). However. Ide Jongl (1 19961) showed 
that this degeneracy c an be partially broken by adding NIR ph otometry to the optical colors, which was 
recently supported bv lWu et al. I (l2005b . ICardiel et al. I (l2003b found that inclusion of an infrared pass- 
band can improve the predictive power of the stellar population diagnostics by ~30 times compared to 
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Table 1 Ultraviolet, optical broad-band, and NIR 2MASS photometry of Gl. 



Filter 


Magnitude (uncertainty) 


Reference 


FUV 


18.972 (0.031) 


Revetal. (2007) 


NUV 


1 8 014 ro 01 ?~l 


u 


14.85 (0.08) 


vandenBerah (1969) 


B 


14.584 (0.013) 


Reed. Harris & Harris (1994) 


V 


13.750 (0.007) 




R 


13.191 (0.010) 




I 


12.684 (0.060) 


Barmbv & Huchra (2001) 


J 


11.858 (0.054) 


Galletietal. (2004) 


H 


11.127 (0.054) 






11.016 (0.054) 





using optical photometry al one. Since NIR photometry is les s sensitive to interstellar extinction than the 
classical optical passbands, iKissler-Patig et alH (l2002l) and IPuzia et al~l (l2002a[) also suggested that it 
provides useful complementary information that can help to disentangle the age-metallicity degeneracy 
(also see Galleti et al. 2004). 

The M31 field was observed as part of the Nearby Galaxies Survey (NGS) by the Galaxy Evolution 
Explorer (GALEX) in two ultraviolet bands (see for details from Rey et al 2005, 2007). lRey et al. I (l2007l) 
pubUshed photometric data for 485 and 273 M31 GCs in the GALEX NUV and FUV bands, respec- 
tively. Gl was detected in these two ultraviolet bands. The GALEX photometric system is calibrated to 
match the spectrophotometric AB system. 

Ivan den Bergh I (1 19691) determined photo-electric photometry f or 45 M31 GCs, including Gl, in the 
UBV bands. Using CCD imaging from the KPNO 0.9m telescope. fReed. Harris & Harris I (Il994 pub- 
lished integrated BVR magnitudes and color indices for 4 1 GCs and GC candidates, including G 1 , in the 
outer halo of M3 1 . We compared the photometry of Gl in the B and V bands bet ween these two studies; 
the results match closely. In this paper, we adopt the CCD BVR photometry of iReed. Harris & Harris I 
(fl994i) . and the photographic ?7-ban d photometry of van den Berg h (19691). with a photometric uncer- 
tainty of 0.08 mag as suggested by iGalleti et al. I (12004 ). Based on HST images. Barmbv & Huchral 
(l200lh detected and published photometry for 1 14 GC candidates associa ted with M31, includin g 32 
new objects. Their F-band p hotometry is in good agreement with that of Ivan den Bergh! d 19691 ) and 
Reed. Harris & Harris ( 19941) . although they do not provide their photometric uncertainties. However, 



Barmbv & Huchra I(l200lh compared their HST photometry with the ground-based measurements com- 



piled by Barmbv et al. ( 2000l) . and found that the median offset in / is 0.06 ± 0.04 mag. Therefore, 
we adopt 0.06 mag as the photometric uncert ainty in the / band. Using the Two Micron All Sky 
Survey (2MASS) database. lOalleti et al. | (|2004 iden t ified 6 93 known and candidate GCs in M31, and 
Usted their 2MASS JHKs magnitudes. iGalleti etaTI (|2004 transformed all 2MAS S magnitudes to th e 
CIT photometric system (Elias et al. 1982, 1983) using the color transformations in lCarpenterl (1200 ll) . 
However, we need the original 2MASS JHKs magnitudes to compare our observa tional SEDs with the 
SSP models, so we reversed this transformation using the same procedures. Since IGalleti et al.1 ( |2004|) 
do not provi de the photo metric uncertainties in JHKs, we obtained these by comparing the magnitudes 
with fig. 2 of lCar penter et al. (2001), where the observed photometric rms uncertainties in the time se- 
ries are shown as a function of magnitude, for stars brighter than the observational completeness limits. 
In fact, th e photo metric uncertainties adopted do not affect our results significantly, as we showed in 
iFan et al."1 (l2006h (see their section 4.3 for details). The full set of ultraviolet, optical broad-band, and 
2MASS NIR ph otometry of Gl is hsted in Table 1. The UBV HI and 2MASS magnitudes are given in 
the Vega svstem ISchneider et al. I (Il977h . For convenience, we converted all observational magnitudes 
to the AB system, following the procedures recommended in BC03. 
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2.3 Reddening and metallicity 

To obtain the intrinsic SED of Gl, the photometric data must first be dereddened. Since Gl is located 
in the halo of M31, i.e., far away from the galaxy's disk, it is (for all practical purposes) only affected 
by Galactic foreground extinctio n. In fact, some auth ors have demonstrated that Gl is affected by a 



negligible amount of reddening. Mevlan et al. I (1200 1[) used HSTAVFPC2 observations in the F555W 
and F814W filters, and applied ISaraiedini I (ll994 's method to simultaneously determine the cluster's 
reddening and metallicity; they obtained a redd ening of E(V — I) = .05 ±0.02 mag toward Gl, which 
is less than the Galactic foreground extinction. Ivan denBerghl ( ll969l) studied the reddening in the halo 
of M31 by comparing the colors of clusters with the same Une-strength index in the Galaxy and in M31, 
and obtained a mean reddening of E{B ~ V) — 0.08 ± 0.02 mag for the clusters in the halo of M31. 
iBarmbv et al. I ( l2000l) determined the reddening for each individual cluster using correlations between 
optical and infrared colors and metallicity, and by defining various 'reddening-free' parameters based 
on their large database of multi-color photometry. For Gl, Barmby et al. (2000, also P. Barmby, priv. 
comm.) obtained E( B — V) = 0.09 ± 0.02 mag. In this paper, we adopt the reddening value from 
iBarmbv et al. I (l2000l) . The val ues for the extinction coefficient, R\, were obtained by interpolating the 
interstellar extinction curve o f lCardelUet ari(ll989h . 

Cluster SEDs are determined by the combination of their ages and metallicities, which is often 
referred to as the age-metallicity degeneracy. Therefore, the age of a cluster can only be constrained 
accurately if the metallicity is known with confidenc e , from independent determinations. There exist 
two metallicity determinations for Gl. iHuchra et al. I ( 199lb derived metallicities for 150 M31 GCs, 



including Gl, using the strengths of six absorption f eatures in the cluster s' integrated spectra. The re- 
sulting metallicity of Gl is [Fc/H] = -1.08 ± O.OQ. lMevlan et al. I ( l200ll ) used HSTAVFPC2 photom- 
etr y to construct deep CM Ds for Gl, combined with the shape of the red-gian t branch as ca l ibrate d 
bv ISaraiedini et al. ] (12000"), to derive the mean metallicity of Gl on the scale of lZinn & West I (11984 . 
[Fe/H] = -0.95 ± 0.09. In this paper, we adopt [Fe/H] = -1.08 ± 0.09 for Gl. 

3 THE STELLAR POPULATION OF Gl 

3.1 Stellar populations and synthetic photometry 

In this section, we compare the SED of Gl with theoretical stellar population synthesis models. We start 
by us ing the BC03 SSP models, which have been upgraded from the earlier iBruzual & Chariot I d 19931 
1 19961) versions, and now provide the evolution of the spectra and photometric properties for a wider 
range of stellar metallicities. BC03 provide 26 SSP models (both of high and low sp ectral resolution) 
using the Padova-1994 evolutionary tracks, half of which were computed based on the lSalpeterl (Il955h 
IMF with lower and upper-mass cut-offs of mr, = 0.1 Mq and rmj = 100 Mq, respectively. The other 
thirteen were computed using the lChabrierl(l2003h IMF with the same mass cut-offs. In addition, BC03 
provide 26 SSP models using the Padova-2000 evolutionary tracks. In this paper, we will use all of these 
SSP models to determine the most appropriate age and mass for Gl. These SSP models contain 221 
spectra describing the spectral evolution of SSPs from 1.0 x 10'' yr to 20 Gyr. The evolving spectra 
include the contribution of the stellar component at wavelengths from 9lA to 160/im. 

Since our observational data are integrated luminosities through a given set of filters, we convolved 
the theoretical SSP SEDs of BC03 with the FUV and NUV, broad-band UBVRI, and 2MASS JHK^ 
filter response curves to obtain synthetic ultraviolet, optical, and NIR photometry for comparison. The 
synthetic magnitude in the AB magnitude system for the itli filter can be computed as 

mx,=~2.5\og -'\. y - 48.60, (1) 

where F\ is the theoretical SED and Lpi the response curve of the zth filter F\ varies as a function of 
age and metallicity. 
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3.2 Fit results 

We use a minimization approach to examine which SSP models are most compatible with the ob- 
served SEDs, following 

X^ = y2-^ 2 (2) 

where TO™°'^(t) is the integrated magnitude in the zth filter of a theoretical SSP at age t, m^^^'^ represents 
the intrinsic integrated magnitude in the same filter, and 

Here, cr^j^g ^ is the observational uncertainty, and crf-^^^^ j is the uncertainty associated with the model 
itself, for the ith filter. ICharlot et al."1(ll996l) estimated the uncertainty associated with the term cr^'^joj j by 
comparing the co lors obtained from different stellar evolutionary tracks and spectral libraries. Following 
IWu et al. 1(120051) . we adopt cr^ = 0.05. 



The BC03 SSP models based on the Padova-1994 evolutionary tracks include six initial 
metallicities, Z = 0.0001,0.0004,0.004,0.008,0.02(^0), and 0.05, corresponding to [Fc/H] = 
-2.25, -1.65, -0.64, -0.33, +0.09, and +0.56. However, the BC03 SSP models based on the Padova- 
2000 evolutionary tracks include six partially different initial metallicities, Z ~ 0.0004, 0.001, 0.004, 
0.008, 0.019 (Zq), and 0.03, i.e., [Fc/H] = -1.65, -1.25, -0.64, -0.33, +0.07, and +0.29. Spectra 
for other metallicities can, in principle, be obtained by interpolation of the appropriate spectra for any 
of these metallicities, although this is not necessarily advisable or straightforward (Frayn & Gilmore 
2002). Instead, we adopt the most appropriate model metallicity for the analysis performed in this pa- 
per. Since we have good estimates of the metallicity and reddening values of Gl (see §2.3), the cluster 
age is therefore the sole parameter to be estimated (for a given IMF and extinction law, both of which 
we assume to be universal). 

None of the SSP models fit the photometric data point in the GALEX FUV band as well as the 
other nine data points. (We checked that the image of Gl in the FUV band is not affected by instrumen- 
tal problems.) Given that Gl contains an old stellar population, the most likely physical explanation for 
this FUV excess compared to the 'standard' BC03 SSP models is the presence of a significant number 
of FUV-bright, hot, extreme horizontal-branch (EHB) stars giving rise to the well-known 'ultraviolet 
upturn' below A ~ 2000A (see, e.g., the review of O'Connell 1999, and references therein; see also 
Landsman et al. 1998; Sohn et al. 2006). (Alternative species, such as AGB-manque stars or blue strag- 
glers are expected to be fewer in number in any 'normal' stellar population.) 

Since 'standard' SSP models do not contain EHB populations, we are forced to deselect the photo- 
metric FUV data point when applying our fitting routines. In Fig. 1, we show the intrinsic SED of Gl 
and the integrated SEDs of the best-fitting models. The dereddened data are shown as the symbols with 
error bars (vertical errors for photometric uncertainties and horizontal error bars for the approximate 
wavelength coverage of each filter); open circles represent the calculated magnitudes of the model SED 
for each filter, obtained by convolving the theoretical SSP SEDs with the appropriate filter response 
curves. The best reduced-x^ values and ages are listed in Table 2. The mass of Gl, also listed in Table 
2, can be estimated by comparing the measured luminosity in the V band with the theoretical mass-to- 
light (M/ L) ratios. These A// L ratios are a function of the cluster age and metallicity. The mass-to-light 
ratios of Gl, calculated based on the metallicity adopted and the age obtained in this paper, are listed 
in Table 2 for the BC03 SSP models. Based on its present luminosity, V = 13.750 ± 0.007 mag, 
and extinction, E{B — V) ~ 0.09 m ag, the cluster's visual magnitude corrected for the extinction is 



Vo = 13.471 ± 0.007 mag, assuming a lCardelU et al. I (Il989h Galactic reddening law with Ay = 0.279 
mag. 

(We note that the NUV data point is also marginally affected by the onset of the UV upturn, which 
causes a slight mismatch between the observations and the best-fitting theoretical SSP models.) 
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Fig. 1 Best-fitting integrated theoretical BC03 SEDs compared to the intrinsic SED of Gl. 
The photometric measurements are shown as the symbols with error bars (vertical for uncer- 
tainties and horizontal for the approximate wavelength coverage of each filter). Open circles 
represent the calculated magnitudes of the model SED for each filter We did not use the FU V 
photometric data point for the fits (see text). 



3.3 Age and mass 

In the previous section we determined the best-fitting age and mass of Gl based on different theoretical 
SSP models. From Table 2 we conclude that, within the errors, the ages obtained from the different 
BC03 models are internally consistent. The mean age of Gl is 18.23 ± 1.76 Gyr. This is in excellent 
agreement with the only other available (rough) age estimate for the cluster by Meylan et al. (2001), 
who estimated its age to be ~ 15 Gyr However, we note that the age of Gl obtained in this paper is older 
than the current-best estimate of the age of the Universe, of order 13.7 Gyr, as discussed in §3.1. We 
will discuss this problem in §4. 

We conclude that the various mass estimates listed in Table 2 place Gl firmly at the top of the cluster 
mass function in the Local Group. Me ylan et al. (200 It) presented three estimates of the total mass of 
Gl, (i) a King-model mass (King 1966) of 1.5 x IO^Mq, (ii) a virial mass of 0.75 x lO'^ M©; and (iii) 
a mass based on a King-Michie model (as defined by Gunn & Griffi 1979) fitted simultaneously to the 
surface brightness profile and the central velocity dispersio n value, estimated be tween 1.4 x lO^M© and 
1.7 X 1O^M0. Our results are in reasonable agree ment with Meylan et alTI ( 1200 lb . although we are aware 



that the King and King-Michie mass estimates of iMeylan et al. ( 2001 ) are up to a factor of two greater 



than ou r photo metric mass estimates. This is not too surprising in view of the model assumptions made. 
ICohen I (|2006|) recently obtained an optical velocity dispersion for the cluster using the Keck/HIRES 
spectrograph, and derived an aperture-corrected line-of-sight velocity dispersion, ctios = 25.5 ±1.5 km 
s^^ (where we have averaged the values she obtained for the two reddest wavelength ranges analyzed; 
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Table 2 Age and mass estimates of Gl based on the BC03 models. 



Evolutionary Track 


IMF 


Age 


/degree of freedom 


M/Lv 


Mass 






(Gyr) 




{M/Lv)q 


(IO^Mq) 


Padova 1994 


Sabeter (-1955) 


19.68 ±0.75 


3.04 


5.10 


1.06 ±0.07 


Padova 1994 


Chabrier ("2003) 


19.79 ±0.50 


2.69 


3.15 


0.65 ±0.04 


Padova 2000 


Sabeter 0955) 


15.44 ±0.78 


5.36 


4.14 


0.86 ±0.05 


Padova 2000 


Chabrier f2003) 


18.01 ± 2.00 


5.27 


2.79 


0.58 ±0.04 



see also Djorgovski et al. 2002). We recently redetermined a projected half-light radius for Gl of ?'h = 
6.5 ± 0.3 pc (Ma et al. 2007b). Thus, based on these most recent results, the dynamical (virial) mass 
of Gl is A/vir = (7.37 ± 2.15) x IO^A/q. This is in excellent agreement with the photometric mass 
estimates obtained in this paper. In turn, this strongly supports the notion that Gl must have had a close- 
to- 'normal' stellar IMF, in order for it to have survived dissolution due to internal two-body relaxation 
until the present time (see also Ma et al. 2006). In particular, this is driven by the observation that if the 
IMF is too shallow, i.e., if a cluster is significantly depleted in low-mass stars compared to (for instance) 
the solar neighborhood, it will disperse within a few orbital periods around its host galaxy's center, and 
most likely within about a Gyr of its formation (e.g., Gnedin & Ostriker 1997; Goodwin 1997; Smith & 
Gallagher 2001; Mengel et al. 2002; Rose, Kou wenhoven & de Grijs, in prep.). 

From the recent work of iMa et al. I ( |2006|) . the intrinsically most luminous M31 GC, 037-B327, 
has been suggested to be the most massive GC in the Local Group, with a total mass of ^ (3.0 ± 
0.5) X 1O^M0, also determined photometrically and somewhat dep ending on the S SP models used, the 
metallicity and age adopted, and the IMF representation. However, ICohen I ( |2006|) pointed out that the 
photometric mass of this cluster had likely been overestimated due to an incorrect extinction correction. 
Nevertheless, she also confirmed the nature of 037-B327 as one of the most massive GCs in the Local 
Group, with a dynamical mass similar to that of Gl. It is intriguing that these two most massive GC in 
M31 bot h are significan tly more massive than the most massive Galactic GC, uo Cen (2.9 — 5.1) x 
106Af^: lMevian1(l2002h l. 

In fact, the high mass of these clusters raises additional, intriguing questions regarding the nature 
of these objects in general, and of Gl in particular (see also Federici et al. 2007; Ma et al. 2007b, and 
references therein). It has been speculated that these objects may be nucleated dwarf galaxies instead of 
genuine GCs. In support of this notion, we point out that Gieles et al. (2006) suggest that there may be 
a physical upper limit to the mass of a star cluster that is not merely the result of size-of-sample effects. 
This maximum mass depends to some extent on the galactic environment; for their example galaxies, 
M51 and the Antennae system, they find a physical upper limit to the stellar mass of ^ (10^ — 10^)AfQ. 
Our values derived for both the photometric and the virial mass of Gl are well above these suggested 
upper mass limits. This may, therefore, provide an additional (although circumstantial) proverbial nail 
in the coffin of Gl as a normal GC. 

3.4 Luminosity of the hot, extreme horizontal-branch stars 

As discussed in §3.3, none of the BC03 SSP models fit the photometric data point in the GALEX 
FUV band as well as the other nine data points. EHB stars may be responsible for this excess in FUV 
band. In fact, iRichet al.l(ll996h found some bluer horizontal-branch stars extending to {V ~ I) = 0.0 
mag, based on their HSTAVFPC2 observations. In this section, we will calculate the luminosity of the 
EHB stars. We assume that the excess in the GALEX FUV band is solely due to these EHB stars. The 
magnitude differences between the four SSP models employed in this paper and the photometric data 
points are 0.63, 1.00, 0.56, and 0.94 mag, respectively, corresponding (2.44, 1.07, 2.11, and 1.00) x 
IQ^Lq, respectively, with a mean number of EHB stars in Gl of (1.65 ± 0.63) x IO^Lq. 
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3.5 Comparison between Gl and S312 

iBrown et al. I (|2004 showed that a 10 Gyr old population in M31 has a main-sequence turnoff at about 
™F8i4W = 28.8 mag. Their deep observations needed exposures of 39.1 hours in F606W and 45.5 
hours F814W, spanning 120 orbits of HST/ACS imaging observations (Brown et al. 2003, 2004). In 
the future, such deep HST observatio ns will only be ob tained for a very small number of fields (e.g.. 
Rich et al. 2005). As discussed in Sl. lMa et al.1 ( l2007al) constrained the age of the M31 GC S312 by 
comparing multicolor p hotometry and the oretical stellar population synthesis models. It is encouraging 
that the age obtained by iMa et al. I (12007 al) is in good agreement with the previous determination based 
on main-sequence photometry (Brown et al. 2004), i.e., 9.5^q gg Gyr versus lOl^'^ Gyr. S3 12 is one 
of the few extragalactic GCs for which the age can be determined from main-sequence photometry. By 
comparing the ages of S3 12 and Gl determined using the same method, we can conclude that S3 12 is 
younger than the majority of the Galactic GCs at the same metallicity, and Gl is as old as the oldest 
Galactic GCs. In fact, if we try to fit the intrinsic SEDs of Gl by the theoretical BC03 SSP SEDs at an 
age of 10 Gyr, the resulting fit is very poor indeed, particularly in the ultraviolet. Therefore, we conclude 
that the method used in this paper, by which the ages of both S3 12 and Gl have been determined, can 
be used to determine the ages of old stellar populations to satisfactory precision, and in particular to 
distinguish between young and old populations. 

4 SUMMARY AND CONCLUSIONS 

In this paper, we first determined the age and mass of the M31 GC Gl, as well as the realistic uncer- 
tainties associated with these estimates, by comparing its multicolor photometry with theoretical stellar 
population synthesis models. Our multicolor photometric data were obtained from GALEX FUV and 
NUV, broad-band optical U BVRI, and 2MASS JHKs observations, which form an SED covering the 
wavelength range from 2267 to 20,000A. Our results confirm that Gl is one of the oldest and most mas- 
sive GCs in the Local Group - that is, if it is indeed a genuine GC given that its mass is well in excess 
of the physical maximum mass predicted by the models of Gieles et al. (2006). 

The age and mass obtained in this paper are somewhat dependent on the SSP model adopted. It is 
evident that the age of 18.23 ± 1.76 Gyr for Gl based on the BC03 models is greater than the currently 
accepted age of the Universe. However, we must keep in mind that the BC03 SSP models were calculated 
for ages up to 20 Gyr. In fact, ages derived for objects such as GCs and galaxies in excess of that of the 
Universe only mean that these objects are among the oldest objects in the Universe. 

In the context of the BC03 models and their associated age range up to 20 Gyr, our derived age for 
Gl is consistent with the suggestion by the WMAP team that the oldest GCs may have formed within 
the first 1.7 Gyr after the Big Bang (see §1). 

The integrated FUV flux depends mainly on the fractional number of horizontal-branch (HB) stars 
with temperatures hotter than ^ 10, 000 K, with a modest dependence on their temperature distri- 
bution (see Rey et al. 2007 and references therein). Older GCs produce more of these hot HB stars and 
they are thus more lik ely to produce str onger FUV fluxes at a given metallicity (see Rey et al. 2007 
and references therein) iLee et al."1 (l2003b showed that the addition of FUV photometry to optical data 
can discriminate cleanly among young (<1 Gyr), intermediate-age (3-5 Gyr), and old (> 14 Gyr) GCs. 
Young and very old GCs exhibit a significant FUV-to-optical spectral continuum slope , but intermediate- 
age clusters are relatively faint in the FUV (see fig. 2 of Lee et al. 2003). Figure 6 of iRey et al. 1 (l2007h 
implies that the age of Gl may be similar to that of the oldest GCs. 

Overall, we therefore conclude that Gl is indeed among the oldest and most massive building blocks 
of M31, and provides a key limitation to the age of the Universe, although we caution that our results 
also provide circumstantial support to the suggestion that the cluster may not be a genuine GC. 
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